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THE  PHASE  SENSITIVITY  OF  AN  INFINITE  LENGTH 
OPTICAL  FIBER  SUBJECTED  TO  A  FORCING  FUNCTION 
AT  A  DEFINITE  FREQUENCY  AND  WAVENUMBER 

1.  INTRODUCTION 

The  phase  sensitivity  of  an  optical  fiber  is  proportional  to  the  longitudinal  and  radial 
strain  in  the  fiber  (Lagakos  and  Bucaro,  1981).  These  strains  are  related  to  the  displacement 
field  of  the  fiber  by  partial  derivatives  of  the  displacements  (Timoshenko  and  Goodier,  1934), 
and  can  be  changed  by  varying  the  external  load  on  the  fiber.  Thus,  the  propagation  of  light 
through  the  fiber  can  be  modified  by  the  application  of  a  mechanical  force  to  the  fiber  exterior. 
In  general,  the  optical  fiber  can  be  modeled  as  an  elastic  solid.  The  displacement  field  of  elastic 
rods  has  been  studied  by  numerous  researchers  because  of  their  use  in  many  mechanical 
designs  (Holden,  1951;  Kaul  and  McCoy,  1964;  Meeker  and  Meitzler,  1964;  Graff,  1975). 
There  has  not,  however,  been  a  thorough  investigation  of  the  phase  sensitivity  of  optical  fibers 
subjected  to  a  forced  response  at  a  definite  frequency  and  wavenumber. 

In  this  report,  the  effects  of  an  infinite  length,  axisymmetric,  isotropic  rod  subjected  to 
an  external  load  at  a  definite  wavenumber  and  frequency  are  presented.  The  resulting 
mechanical  stresses  are  used  in  fiber  optic  equations  to  determine  the  phase  sensitivity  and  the 
change  of  refractive  index  of  the  fiber.  The  derivation  begins  with  the  partial  differential 
equations  of  motion  of  an  isotropic,  elastic,  slightly  damped  solid  in  cylindrical  coordinates. 
The  displacements  are  written  as  the  sum  of  a  dilatational  component  and  an  equivoluminal 
vector  component.  This  approach  allows  the  equations  to  be  separated  into  a  dilatational  wave 
equation  and  three  distortional  (shear)  wave  equations.  General  solutions  to  these  four  wave 
equations  are  next  determined.  These  general  solutions  are  inserted  at  the  diameter  of  the  rod 
into  the  stress- strain  relationships,  which  are  equated  to  the  externally  applied  loading  function 
on  the  cylinder.  Based  on  these  relationships,  a  two-by-two  system  of  linear  equations  is 
developed.  The  solution  to  the  linear  equations  yields  the  displacements  of  the  rod.  The 
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strains  are  determined  from  these  displacements  and  then  inserted  into  equations  that  yield  the 
optical  phase  sensitivity  and  the  change  of  refractive  index  of  the  fiber. 


2.  SYSTEM  MODEL  AND  CLOSED-FORM  SOLUTION 

The  system  model  is  a  cylindrical,  linear,  isotropic  medium  whose  motion  is  governed 
by  the  equation  (Timoshenko  and  Goodier,  1934) 


+  (A  +  /t)  VV  •  u  =  p 


(1) 


where  p  is  the  density;  A  and  p  are  the  Lamd  constants;  t  is  time;  •  denotes  a  vector  dot 
product;  u  is  the  cylindrical  coordinate  displacement  vector  expressed  as 


Ur(r,d,z,t) 
u  =  \uQ(r,d,z,t) 


(2) 


with  subscript  r  denoting  the  radial  direction,  d  denoting  the  angular  direction,  and  z  denoting 
the  axial  direction;  V  is  the  gradient  vector  differential  operator  written  in  cylindrical 
coordinates  as  (Potter,  1978) 


^  d  .  Id.  d  . 


(3) 


with  ij.  denoting  the  unit  vector  in  the  r  direction,  Iq  denoting  the  unit  vector  in  the  6  direction, 
and  /  denoting  the  unit  vector  in  the  z  direction;  is  the  three-dimensional  Laplace  operator 


operating  on  vector  u  as 


y2„  _ifL_ 

V  Ur  2 


2  dug  Y 


r2  r2  dd  F 


Ir  + 


with  operating  on  scalar  u  as 


Y72„  -  V  ,  ^  ^  {  ^r,e,z  1  . 


(4) 


(5) 


and  the  term  V  •  u  is  called  the  divergence  and  is  equal  to 
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_  dur  1  dufi 
V*u  =  ^  +  -^  +  -T-^  +  —  . 
3r  r  99  dz  r 


dUj  Ur 

dz 

The  coordinate  system  of  the  rod  is  shown  in  figure  1. 


Figure  1.  Cylindrical  Rod 


(6) 


The  displacement  vector  u  is  written  as 

u  =  V^  +  VxH,  (7) 

where  ^  is  a  dilatational  scalar  potential,  x  denotes  a  vector  cross  product,  and  H  is  an 
equivoluminal  vector  potential  expressed  as 


H  = 


Hr{r,d,z,t) 

HQir,e,z,t) 


(8) 


Expanding  equation  (7)  and  breaking  the  displacement  vector  into  its  individual  terms  yields 
^  dr  r  dd  dz  ' 
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and 


and 


_  1  d(p  dH,  ^ 
rdd^  dz  dr 


_d(l>HQ  dHe  1  dHr 

^  dz  r  dr  r  dd 


Equation  (7)  is  next  inserted  into  equation  (1),  which  results  in 


"  dt^ 


(10) 


(11) 


(12) 


(13) 


The  constants  and  are  the  complex  dilatational  and  shear  wave  speeds,  respectively,  and 
are  determined  by 


Cd=. 


I A  -("  2ii 


(14) 


and 


Ci  =  I—  • 


(15) 


The  relationship  of  the  Lamd  constants  to  the  compressional  and  shear  moduli  is  shown  as 
Ev 


A  =  - 


(l  +  u)(l-2u) 


(16) 


and 


/i  =  G  = 


E 

2(1 +  t>)  ’ 


(17) 


where  E  is  the  complex  compressional  modulus  (N/m^),  G  is  the  complex  shear  modulus 
(N/ra^),  and  u  is  the  Poisson's  ratio  of  the  material  (dimensionless). 

The  conditions  of  infinite  length,  axisymmetric  response  (n  =  0),  and  steady- state 
response  are  now  imposed,  allowing  the  scalar  and  vector  potential  to  be  written  as 
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(p  =  g{r)cosind)e^^"e‘^  =  g{r)e‘^e''^  , 

(18) 

Hr  =hf.{r)sin{nd)e‘'^^e^^  =0  , 

(19) 

Hg  =hg(r)cos(nd)e^^e^^  =hg(r)e'’^^e‘^  , 

(20) 

Hz  =  h7(r)sm(nd)e^^e^^  =0  , 

(21) 

where  k  is  the  wavenumber  of  excitation,  (O  is  the  frequency  of  excitation,  and  i  is  the  square 
root  of  -1.  Note  that  for  axisymmetric  response,  the  equations  of  motion  are  dependent  only 
on  the  scalar  potential  (p  and  angular  contribution  Hq  of  the  vector  potential.  Additionally, 
because  Hj.  =  0,  =  0,  and  Hg  and  (p  are  not  functions  of  6,  equation  (10)  becomes 


.,-§-0.  (22, 

where  (  )  denotes  any  function.  Inserting  equations  (18)-(21)  into  equations  (12)  and  (13) 
yields  the  following  four  wave  equations: 


d^gjr)  1  dg(r) 

dr^  r  dr 


ycl 


8{r)  =  Q  , 


(23) 


d\{,r)  ,  1  dhr{r)  , 
dr^  r  dr 


(O 


1 


A 


V 


2  ^  2 
r 


hr{r)  =  Q  , 


d%{r)  .  1  dhgir)  . 
dr^  r  dr 


VC? 


hg{r)  =  0  , 


and 


(24) 


(25) 


d^h^ir)  .  1  dh^ir)  . 

^  1  ^  "T 


dr^ 


r  dr 


V 


\ 


hAr)  =  Q  . 


(26) 


J 


The  solution  to  equations  (23)  and  (25)  is  now  found  with  a  Bessel  function.  No 
solution  is  found  to  equations  (24)  and  (26)  because  they  do  not  contribute  to  the  axisymmetric 
response.  The  solution  to  equation  (23)  is 
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g{r)  =  CiJQ{ar)  ,  (27) 

where  Jq  is  a  complex,  zero-order,  first-kind  standard  Bessel  function;  Cj  is  a  complex 

constant  (determined  below);  and 


a  = 


(28) 


The  solution  to  equation  (23)  is 

he{r)  =  C2hm  .  (29) 

where  is  a  complex,  first-order,  first-kind  standard  Bessel  function;  C2  is  a  complex 
constant  (determined  below);  and 


(0^ 


(30) 


The  numerical  evaluation  of  the  complex-valued  Bessel  functions  is  described  in  the  appendix. 

The  displacements  and  external  forces  are  now  equated  by  use  of  the  stress-strain 
constitutive  equations  on  the  forced  surfaces  of  the  rod.  The  normal  stress,  strain,  and  radial 
forces  acting  on  the  cylinder  are  related  by 

CTrr (a, 0, z, 0  =  (2-  +  2 ^)err (a, 0, z, t)  -f-  ?^Qe («’ 0  +  ^zz t)  =  p{a,e,z,t)  ,  (31) 


where  <Trr(a,0,z,t)  is  the  normal  radial  stress,  £rr(a,0,z,t)  is  the  normal  radial  strain, 
£QQ(a,d,z,t)  is  the  normal  circumferential  strain,  £^{a,d,z,t)  is  the  normal  longitudinal  strain, 
p{a,d,z,t)  is  the  external  pressure  on  the  cylinder  in  the  radial  direction,  and  a  denotes  the 
outer  radius.  The  shear  stress,  strain,  and  longitudinal  forces  are  related  using 

CTrz  (a,  d,  z,  t)  =  2p£rz  (a,  9,  z,  t)  =  f{a,  9,  z,  t)  ,  (32) 

where  (7rz{a,9,z,t)  is  the  shear  stress,  £rz{a,9,z,t)  is  the  shear  strain,  and  fia,9,z,t)  is  the 
external  shear  stress  on  the  rod  in  the  longitudinal  direction. 

The  strains  are  related  to  the  displacements  in  an  axisymmetric  solid  by  (Timoshenko 
and  Goodier,  1934) 
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^Ur 

(33) 

ew=J  . 

(34) 

duj 

(35) 

and 

1  ( du^  duj-  \ 

(36) 

The  relationship  between  the  displacements  (and  the  derivatives  of  the  displacements)  and  the 
potential  functions  g  and  hg  is  found  by  combining  equations  (9),  (11),  (18),  and  (20)  to 

produce 


(37) 

(38) 


duj-  _ 
dr 

du^ 
dz 

du^ 
dr 


(  a1 


d  g(r)  -j-dh^jr) 


dr^ 


dr 


ikZeiox 


r  dr  J 


.,  dg(r)  ,  1  dhgir)  hgjr)  .  d^hgjr) 


ik-^  +  - 
dr  r  dr 


dr^ 


ikZeiox 


(39) 

(40) 

(41) 


and 


^  =  r ikMll  +  k\(r)ykz^^^  .  (42) 

dz  \  dr  J 

Combining  equations  (31),  (33),  (34),  (35),  (37),  (39),  and  (40)  yields  the  normal 
stress  in  terms  of  the  potential  functions  g  and  hg  air  =  a  as 

(A  +  2n)^^  +  -  M^gia)  -  =  P  ,  (43) 

dr^  a  dr  dr^ 
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where  P  is  the  magnitude  of  the  normal  force  acting  on  the  exterior  of  the  cylinder.  Combining 
equations  (32),  (36),  (41),  and  (42)  yields  the  shear  stress  in  terms  of  the  potential  functions  g 
and  at  r  =  a.  The  result  is 


lidhQia)  ,  „d^hQ{a)  _  ^ 

j  ^  ”  j  2  “  ^  ’ 

a  dr  dr 


(44) 


where  F  is  the  magnitude  of  the  externally  applied  shear  stress  acting  on  the  exterior  of  the  rod. 
Implicit  in  equations  (43)  and  (44)  is  the  assumption  that  the  external  loads  on  the  cylinder  are 
occurring  at  a  definite  frequency  and  wavenumber. 

Inserting  equations  (27)  and  (29)  into  equations  (43)  and  (44),  applying  the  recurrence 
relationships  of  the  first-kind  standard  Bessel  function,  and  then  rewriting  as  a  two-by-two 
system  of  linear  equations  results  in 


an  012 

L^21  ^22 J 


Q; 

,C2J 


If( 


(45) 


where  the  matrix  coefficients  a^^j  are  given  as 


Oil  =  j^(-A  - 2^)0^  -  j7o(ow)  + 

f  2jj.ik' 


\  a 


and 


ai2=(-2/iWoW  +  l' 

a2i  ■=  {-2iiika)Ji{aa)  , 


hm , 


(46) 

(47) 

(48) 

(49) 


The  radial  displacement  of  the  rod,  now  found  using  equations  (27),  (29),  and  (37)  and 
the  constants  Cj  and  C2  from  equation  (45),  is  given  as 

ur{r,z,t)  =  [-Cia7i(ar)  -  C2ifc/i()3r)>‘V“'  , 
which  can  also  be  written  as 


(50) 


Ur(r,z,t)  =  Ur(r)e‘^e‘^ 


(51) 
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The  longitudinal  displacement  of  the  cylinder  is  found  using  equations  (27),  (29),  and  (38), 
and  is  written  as 

U,(r,z,t)  =  [CiiUQiar)  +  (52) 

or 

u,{r,z,t)  =  U,{r)e^’^V^  .  (53) 

The  strains  are  evaluated  at  r  =  0  because  light  propagation  is  confined  primarily  to  the  middle 
of  the  fiber.  The  radial  strain  of  the  rod  from  equation  (33)  is 

=  .  -  (54) 

L  ^  ^ 

The  longitudinal  strain,  found  from  equation  (35),  is 

(0,  z,  0  =  [-Cik^  +  C2  (ikp)]  .  (55) 

The  strains  are  related  to  the  optical  phase  sensitivity  by  the  equation 

^  =  %-^[£rrWl  +  f!2)  +  %/i2]  .  (56) 

(p  2 

where  hq  is  the  refractive  index  of  the  fiber,  and  ^11  and  Pj2  Pockels'  coefficients. 

Inserting  the  values  of  /iq  ~  T46,  =  0.126,  and  P22  =  0-270  changes  equation  (56)  to 

^  =  0.112£^  -  QA22£rr  •  (57) 

The  strains  are  related  to  the  change  of  refractive  index  by  the  equation 

_„3 

dn  =  -^[2£rr{Pi-PA)  +  ezz{Pi-^PA)]  -  (58) 

where  P\  -Pn  -  (^-  126  and  P4  =  (P11-P12)  /  2  =  -0.0718.  Inserting  these  values  changes 
equation  (58)  to 

= -0.616  e^r- 0.420  £^2  .  (59) 
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The  above  equations  show  that  given  the  physical  properties  of  the  fiber  and  an  external  forcing 
function  at  a  definite  wavenumber  and  frequency,  the  optical  phase  sensitivity  and  the  change 
of  refractive  index  can  be  calculated  in  a  closed-form  solution. 

3.  FINITE  ELEMENT  SOLUTION 

The  closed-form  solutions  (equations  (57)  and  (59))  are  compared  to  results  obtained 
with  finite  element  analysis  for  the  specific  case  of  an  external  normal  force  applied  at  zero 
wavenumber.  Finite  element  analysis  is  a  discretized  modeling  technique  that  breaks  down  the 
structure  into  a  number  of  subdomains  (elements)  (Cook,  1974;  Zienkiewicz,  1983). 
Constitutive  equations  are  applied  to  each  of  these  elements  to  develop  equations  of  motion. 
The  individual  equations  of  motion  are  then  assembled  to  form  a  global  equation: 

M  +  C  +  Kg(t)  =  f(x,0  ,  (60) 

dt 

where  M  is  the  system  mass  matrix,  C  is  the  system  damping  matrix,  K  is  the  system  stiffness 
matrix,  g(0  is  the  generalized  displacement  vector,  and  f(x,r)  is  the  external  forcing  function 
vector.  The  finite  element  results  displayed  here  were  obtained  using  ANSYS,  which  is  a 
general  purpose  finite  element  computer  program.  Finite  element  analysis  has  the  advantage  of 
handling  material  and  property  discontinuities  that  are  not  easily  resolved  with  other  numerical 
techniques.  It  is  routinely  used  for  comparisons  with  other  experiments  and  theories  to  verify 
results.  It  is  important  to  note  that,  although  the  finite  element  analysis  is  an  approximate 
numerical  technique,  it  is  extremely  different  from  the  analytical  solution  presented  in  section  2. 

The  finite  element  grid  of  the  fiber  is  shown  in  figure  2.  Because  the  elements  in  this 
analysis  are  axisymmetric,  the  torsional  motion  of  the  rod  is  zero.  The  problem  was  discretized 
with  10  elements  in  the  radial  direction  and  was  loaded  with  a  normal  forcing  function 
(harmonic  in  time)  on  the  exterior  of  the  cylinder.  The  forcing  function  had  no  spatial  variation 
and  therefore  was  equivalent  to  zero  wavenumber  excitation.  The  outputs  of  the  finite  element 
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program  are  longitudinal  and  radial  stress;  these  values  were  inserted  into  equations  (57)  and 
(59)  for  comparison  of  the  finite  element  method  to  the  closed-form  solution. 


Figure  2.  Finite  Element  Grid 


4.  ANALYSIS 

An  analysis  of  an  optical  fiber  with  the  above  mathematical  methods  was  performed. 
The  following  mechanical  parameters  were  used:  a  =  62.5  x  10’^  m,E-  7.2  x  10^^  (1  -h 
0.010  N/m^,  u  =  0.17  (dimensionless),  and  p  =  2200  kg/m^.  The  test  case  of/=  50  Hz,  k  = 
0,  F  =  0,  and  F  =  1  N/m^  is  used  for  model  verification.  These  results  are  listed  in  table  1. 
The  experimental  results  in  the  table  are  from  a  previous  study  (Davis  et  al.,  1982).  Note  that 
all  three  methods  show  results  within  0.16  dB  of  each  other. 


Table  1.  Comparison  of  Results  at  Zero  Wavenumber 


Method 

1  Ap/  pF  1  (m^/N) 

|A^/pFI(dBrefm2/N) 

Closed-Form  Model 

4.54  X  10-12 

-226.86 

Finite  Element  Model 

4.58  X  10-12 

-226.78 

Experiment 

4.50  X  10-12 

-226.94 
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The  wavenumber  analysis  of  the  model  was  conducted  at  two  frequencies:  50  and  500 
Hz.  This  analysis  was  accomplished  by  first  setting  F  =  0  and  P  =  1  N/m^  to  determine  the 
optical  response  of  the  fiber  to  a  normal  pressure  excitation  of  unity  and  by  next  setting  F  =  1 
N/m^  and  F  =  0  to  calculate  the  optical  response  of  the  fiber  to  a  shear  force  excitation  of  unity. 
In  figures  3-18,  the  short  dashed  line  represents  the  dilatational  wavenumber  {k^  =0) I c^i)  and 
the  long  dashed  line  represents  the  shear  wavenumber  (kg=0)l  c^) .  The  exponential  functions 

in  time  and  space  that  are  associated  with  the  closed-form  solution  are  suppressed. 

Figure  3  is  a  plot  of  the  transfer  function  of  longitudinal  strain  divided  by  normal 
pressure  versus  wavenumber  for/  =  50  Hz.  Figure  4  is  a  plot  of  the  transfer  function  of  radial 
strain  divided  by  normal  pressure  versus  wavenumber  for/=  50  Hz.  Figure  5  is  a  plot  of  the 
transfer  function  of  optical  phase  sensitivity  divided  by  normal  pressure  versus  wavenumber 
for/=  50  Hz.  Figure  6  is  a  plot  of  the  transfer  function  of  change  of  refractive  index  divided 
by  normal  pressure  versus  wavenumber  for  /  =  50  Hz.  Figure  7  is  a  plot  of  the  transfer 
function  of  longitudinal  strain  divided  by  shear  stress  versus  wavenumber  for/=  50  Hz. 
Figure  8  is  a  plot  of  the  transfer  function  of  radial  strain  divided  by  shear  stress  versus 
wavenumber  for  /  =  50  Hz.  Figure  9  is  a  plot  of  the  transfer  function  of  optical  phase 
sensitivity  divided  by  shear  stress  versus  wavenumber  for/  =  50  Hz.  Figure  10  is  a  plot  of  the 
transfer  function  of  change  of  refractive  index  divided  by  shear  stress  versus  wavenumber  for/ 
=  50  Hz. 

Figure  11  is  a  plot  of  the  transfer  function  of  longitudinal  strain  divided  by  normal 
pressure  versus  wavenumber  for/=  500  Hz.  Figure  12  is  a  plot  of  the  transfer  function  of 
radial  strain  divided  by  normal  pressure  versus  wavenumber  for/=  500  Hz.  Figure  13  is  a 
plot  of  the  transfer  function  of  optical  phase  sensitivity  divided  by  normal  pressure  versus 
wavenumber  for/=  500  Hz.  Figure  14  is  a  plot  of  the  transfer  function  of  change  of  refractive 
index  divided  by  normal  pressure  versus  wavenumber  for/  =  500  Hz.  Figure  15  is  a  plot  of 
the  transfer  function  of  longitudinal  strain  divided  by  shear  stress  versus  wavenumber  for/= 
500  Hz.  Figure  16  is  a  plot  of  the  transfer  function  of  radial  strain  divided  by  shear  stress 
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versus  wavenumber  for  /  =  500  Hz.  Figure  17  is  a  plot  of  the  transfer  function  of  optical 
phase  sensitivity  divided  by  shear  stress  versus  wavenumber  for  /=  500  Hz.  Figure  18  is  a 
plot  of  the  transfer  function  of  change  of  refractive  index  divided  by  shear  stress  versus 
wavenumber  for/=  500  Hz. 

Comparing  the  fiber  response  to  normal  excitation  (figures  3-6  and  figures  11-14)  and 
the  fiber  response  to  shear  excitation  (figures  7-10  and  figures  15-18)  shows  that  the  fiber  is 
120  dB  more  sensitive  to  shear  excitation  than  to  normal  excitation  at  50  Hz  and  is  100  dB 
more  sensitive  to  shear  excitation  than  to  normal  excitation  at  500  Hz.  The  response  to  normal 
excitation  has  relatively  little  change  from  50  to  500  Hz.  This  implies  that  the  time  or 
frequency  dynamics  associated  with  a  normal  load  on  the  fiber  (and  a  corresponding  response 
at  r  =  0)  are  extremely  small.  This  was  confirmed  by  analyzing  the  finite  element  model  with 
static  excitation,  which  produced  results  identical  to  a  dynamic  excitation  at  frequency /=  50 
Hz.  The  fiber  response  to  shear  excitation  is  a  decreased  response  of  20  dB  from  50  to  500 
Hz.  In  all  the  figures,  the  maximum  response  is  occurring  very  close  to  the  dilatational 
wavenumber.  Additionally,  the  response  around  the  shear  wavenumber  is  flat,  indicating  that 
there  is  no  shear  wave  propagation  in  the  structure. 

The  above  closed-form  solution  (equations  (l)-(59))  can  be  simplified  for  low- 
wavenumber  regions  with  linear  approximations  to  the  Bessel  functions.  Use  of  the  small 
argument  assumptions  allows  the  Bessel  functions  to  be  written  as 

(61) 

and 


w-f. 

Using  equations  (61)  and  (62), 


(62) 

the  matrix  coefficients  in  equation  (45)  can  be  rewritten  as 


ail 


a\l  =  , 


(63) 

(64) 
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and 


(65) 


a  (66) 

II  2  2 

The  strains  are  evaluated  using  equations  (54)  and  (55),  with  the  constants  Cj  and 
determined  by  equations  (63)-(66).  This  expression  of  the  closed-form  solution  is  slightly 
different  than  that  derived  earlier  because  it  contains  the  low-wavenumber  assumption.  The 
low-wavenumber  solution  was  compared  to  the  exact  closed-form  solution.  The  results  are 
nearly  identical  up  to  a  wavenumber  of  1000  rad/m  and  vary  only  by  0.2  dB  at  a  wavenumber 
of  4000  rad/m  (for  calculations  at  50  Hz).  For  most  applications,  the  low- wavenumber 
assumption  is  adequate  to  describe  the  response  of  the  fiber. 


5.  CONCLUSIONS 

The  optical  phase  sensitivity  of  an  infinite  length  fiber  is  shear  stress  dominated.  When 
the  fiber  is  subjected  to  a  normal  pressure,  the  frequency  domain  dynamics  are  extremely 
small.  When  the  fiber  is  subjected  to  a  shear  stress,  there  is  a  20-dB  drop  in  the  optical  phase 
sensitivity  between  50  and  500  Hz.  The  low- wavenumber  assumption  is  accurate  to  model  the 
response  of  the  fiber  up  to  a  wavenumber  of  at  least  4000  rad/m. 
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Figure  3.  Transfer  Function  of  Longitudinal  Strain  Divided  by  Normal 
Pressure  Versus  Wavenumber  for/=  50  Hz 
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0  0.1  0.2  0.3  0.4 

Wavenumber,  k  (rad/m) 

Figure  4.  Transfer  Function  of  Radial  Strain  Divided  by  Normal  Pressure 
Versus  Wavenumber  for /=  50  Hz 
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Figure  5.  Transfer  Function  of  Optical  Phase  Sensitivity  Divided  by  Normal 
Pressure  Versus  Wavenumber  for/=  50  Hz 
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Figure  7.  Transfer  Function  of  Longitudinal  Strain  Divided  by  Shear  Stress 
Versus  Wavenumber  for/=  50  Hz 
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Figure  8.  Transfer  Function  of  Radial  Strain  Divided  by  Shear  Stress  Versus 

Wavenumber  for/=  50  Hz 
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Figure  9.  Transfer  Function  of  Optical  Phase  Sensitivity  Divided  by  Shear 
Stress  Versus  Wavenumber  for  /=  50  Hz 
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Figure  12.  Transfer  Function  of  Radial  Strain  Divided  by  Normal  Pressure 
Versus  Wavenumber  for/=  500  Hz 
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Figure  14.  Transfer  Function  of  Change  of  Refractive  Index  Divided  by 
Normal  Pressure  Versus  Wavenumber  for  /  =  500  Hz 
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Figure  18.  Transfer  Function  of  Change  of  Refractive  Index  Divided  by  Shear 
Stress  Versus  Wavenumber  for  /=  500  Hz 
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APPENDIX  -  NUMERICAL  EVALUATION  OF  THE  BESSEL  FUNCTION  J 


The  complex  Bessel  function  J  was  evaluated  numerically.  Based  on  the  magnitude  of 
the  argument,  this  process  used  one  of  two  equations.  If  the  magnitude  of  the  argument  of  z 
was  less  than  3.0,  the  equation 


\n  16  /  /  /j  \A: 

=  I'  ’ 


k=0 


k\{k  +  n)\ 


(A-1) 


was  used.  When  the  magnitude  of  the  argument  of  z  was  greater  than  3.0,  an  asymptotic 
expansion  was  used  (where  |argz|  <7r): 


Jniz)  = 


A„(z)cos! 


z  - 


nn 

~2 


Bn{z)sm 


z- 


nK 


(A-2) 


where 


A„(z)  =  l- 


(4n^  -l)(4n^  -9)(4n^  -25)(4n^  -  49)(4n^  -  81)(4n^  - 121) 

6!(8z)^ 


(4n^  -l)(4n^  -9)(4n^  -25)(4n^  -49)(4«^  -81)(4«^  -121)(4n^  -169)(4n^  -225) 
^  8!(8z)^ 


(A-3) 


and 


„  ,  ,  (4n^  - 1)  (4n^  -  l){4n^  -  9)(4n^  -  25) 

3!(8z)’ 


(4n^  -  l)(4n^  -  9)(4n^  -  25)(4ai^  -  49)(4n^  -  81) 
5!(8z)^ 


(4n^  -l)(4n^  -9)(4«^  -25)(4«^  -49)(4n^  -81)(4n^  -121)(4n^  -169) 

7K&? 


(A-4) 


A-1 


INITIAL  DISTRIBUTION  LIST 


Addressee 

Defense  Technical  Information  Center 
Naval  Research  Laboratory  (T.  Dandridge,  N.  Lagakos) 
Naval  Sea  Systems  Command  (Capt.  G.  Kent  (PMS-425), 
D.  Lechner  (PMS-425 1 1)) 


No.  of  Copies 
12 
2 

2 


Office  of  Naval  Research  321  (T.  Goldsberry,  R.  Varley,  K.  Dial,  S.  Littlefield, 


D.  Davison) 

Program  Executive  Office,  USW/ASTO  (Cdr.  J.  Polcari,  W.  Chen,  A.  Hommel, 
J.  Jones,  Lcdr.  M.  Tra week,  R.  Melusky) 

Space  And  Naval  Warfare  Systems  Command  (J.  FeuiUet  (PMW-182)) 


